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Abstract 

In our nonperturbative lattice investigation we study the interface tension of the finite- 
temperature electroweak phase transition. In this analysis the Higgs mass has been chosen 
to be about 35 GeV. At the transition point of a finite volume system, tunnelling between 
the symmetric and the Higgs phase takes place. This phenomenon leads to a splitting of 
the ground state, which can be used to determine the interface tension. The result obtained 
this way agrees with the result of the two-coupling method and with the prediction of the 
perturbative approach. 
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1 Introduction 



At high temperatures the electroweak symmetry is restored. In the early universe a cosmological 
phase transition took place between the high-temperature -"symmetric"- and low-temperature 
-"Higgs"- phases. Since the baryon violating processes are unsuppressed at high temperatures, 
the observed baryon asymmetry of the universe has finally been determined at the electroweak 
phase transition M. 

In recent years quantitative studies of the electroweak phase transition have been carried out 
by means of resummed perturbation theory and lattice Monte Carlo simulations @— 0. 

It has been clarified that in the SU(2)-Higgs model for Higgs masses (m#) below 50 GeV, the 
phase transition is of first order. However, its strength rapidly decreases for increasing Higgs 
boson masses. 

One of the most important quantities that characterize the strength of the phase transition 
is a, the interface tension between the two phases. The course of the cosmological electroweak 
phase transition, thus the nucleation rate and temperature, is substantially influenced by a. 

The determination of o for m# ~ 35 GeV in the SU(2)-Higgs model is the main goal of the 
present letter. We propose and apply some new methods different from those of || . Comparing 
the results we strengthen the confidence in the correctness of the determination of a. 

The plan of the paper is as follows. In Section 2 we give some details of the simulations and 
suggest a new method to determine the transition point. Section 3 contains the determination of 
a based on the finite volume tunnelling phenomena between the two phases. In Section 4 a brief 
description of the two-coupling study of the interface tension with a bootstrap error-analysis is 
presented. Section 5 is devoted to our conclusions. 



2 Simulation 



The lattice action of the SU(2)-Higgs model is given by 



S[U,<p] 



0£(i 

P i 



P l 



+ £ |Tr(^) + A[±Tr(^)-l 



(1) 



where U X;IX denotes the SU(2) gauge link variable and U p i is the path-ordered product of the 
four U x ^ around a plaquette. The symbol ip x stands for the Higgs field, which is also written as 
<fx — Px • ol x , with p x G R + and a x G SU(2). Details of the simulation algorithms can be found 

in g|ig. 



The order parameters under study are p 2 , s\ og , the density of S\ og = S[U,(p] — 3 J2x l°g(Ar)> 
and Lip, the density of = ^xn fTr (fx+pUx^ <Ac)- The latter plays an important role 

in our reweighting [|TTJ, which has been used to obtain data in the vicinity of the simulation 
point Pf. 



In our finite-temperature simulations we have used elongated lattices: V 



Lt x L 2 xy x L z , 



where L t 



2, L xy L z . Keeping the bare parameters /? = 8.0, A = 3.0 - 10 and L t 



2 fixed 



2 



we tuned k to its critical value. These parameters corresponds to a Higgs boson mass of roughly 
35 GeV at T = 0. 

As will be shown in Section 3, the calculation of a by use of the two-coupling method presup- 
poses a precise determination of the critical hopping parameter k c on large lattices. We suggest 
here a new method, which in practice gives the same and very precise k c as the multicanoni- 
cal method |13|]. However, using this method the notoriously large autocorrelation time of the 
multicanonical simulation can be reduced dramatically. 

The method is based on the observation that, for sufficiently large elongated lattices, a flat 
regime appears between the peaks in the probability distribution of an order parameter ||. The 
peaks correspond to pure phase configurations. Between the peaks the system is dominated by 
configurations with two interfaces perpendicular to the long direction z, and two bulk phases. At 
the transition point the free energy densities of these bulk phases are the same, the ratio of the 



occupied volumes is arbitrary, thus the distribution is flat. By use of constrained simulation [12] 
one can enforce that an order parameter of the system stays in a given short interval between 
the peaks. Close to the transition point one can use a reweighting procedure to determine the 
value of the hopping parameter at which the distribution is flat. 

In practice we have chosen si og as the order parameter of the constrained simulation. Since s\ og 
remains unchanged by the overrelaxation algorithms [|l(|, only the heatbath algorithms have to 
be modified. Every proposal that leads to expectation values of sy og outside the selected interval, 
is rejected. The observed acceptance rate of the algorithm is larger than 98%. The lattice size 
under study was 2 x 24 2 x 256. In a series of two short runs without any bounds on s\ og , the 
position of the two peaks, belonging to the pure phases, was estimated. As a consequence the 
interval s\ og G [4.90, 4.95] has been chosen to determine k c . We have performed 18000 measuring 
sweeps at k — 0.12865, and by using reweighting we have obtained results in the vicinity of 
this point. The integrated autocorrelation time turned out to be r int = 21(2) sweeps for s\ og 
and L^^. This value is much smaller than a typical autocorrelation time for a multicanonical 
simulation. The reason is the following. In the case of a multicanonical simulation the system 
must completely go from one phase to the other in order to obtain an independent configuration. 
In our method only the ratio of the bulk phases changes slightly. However, this is enough to tell 
whether the free energy densities are the same or not. 

The transition point is defined by the implicit relation for k c 

(siog)re c = | ( s log,- + Slog,+) > (2) 

with si ogi _ and si og ,+ denoting the bounds on si og . In our case this definition is equivalent to the 
condition that at k c the distribution is flat in the given interval. 

For the determination of the statistical errors on (s\ og ) K we have used the results of a bootstrap 
procedure with 3000 iterations and 30 independent subsamples (see Section 4). The final result 
of this analysis is k c = 0.1286565(7). Within the error bars, it agrees with the result from the 
2 x 16 2 x 256 lattice (see Section 3). Thus, the finite-size effects of k c beyond L xy = 24 are 
expected to be small. The values of 

Slog ^ 

function of k are plotted in figure |1|. 
The above procedure gives the same k c for the upper half and for the lower half of the interval. 
This fact can be considered as a check that the system was in the flat regime. We frequently 
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lattice: 2 x 24 x 24 x 256 
(3 = 8.0 A = 0.0003 




0.1286555 0.1286565 0.1286575 



Figure 1: The order parameter S Jog 3b function of the hopping parameter. The dashed lines 
give the range of the statistical error. From 4.90 < si og < 4.95 the result k c = 0.1286565(7) is 
obtained. 



tested that the configurations contained two bulk phases and two interfaces in between. 

As a non-trivial cross-check of the method we should mention that, for A = 0.0001 (m-H ~ 18 
GeV) the results of || obtained by the multicanonical method have been reproduced with this 
new procedure. 



3 Tunnelling and energy splitting 

The electroweak symmetry is spontaneously broken; however, this phenomenon does not appear 
in lattice simulations, since at finite volume the symmetry is never broken. In this case a small 
energy splitting (Eq) appears between the ground state and the first excited state, which vanishes 
in the infinite volume limit. The appearance of this Eq at finite volume is due to the tunnelling 
effect between the symmetric and the Higgs phase. In the z-direction different domains can 
be observed, with domain walls between them. The associated interface energy determines the 
tunnelling mass (E ) via 

E = Cex P (-L t L% / (r). (3) 
This tunnelling effect has been used to determine, for instance, the interface tension of the 



four-dimensional Ising model in the broken phase [13, [L5|], to calculate the interface tension 
in three-dimensional systems [|16|], to determine the interface tension of the finite-temperature 
confinement-deconfinement phase transition of the pure SU(3) gauge theory [17j], and to deter- 



mine the interface tension of the electroweak phase transition ||18[| . In general, the prefactor C 
of eq. (3) can depend on L xy ; but in our case the leading contribution is L .^-independent. 

According to eq. (0) we have measured the energy splitting on different lattice sizes. The 
number of sweeps and the critical hopping parameters are listed in table |l[ The observed expo- 
nential autocorrelation times are 0(100) in the case of smaller volumes {L xy = 4, 8) and (9(1000) 
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2x12x12x128 


512000 


0.128663(2) 


2x12x12x256 


512000 


0.128665(2) 


2x16x16x128 


448200 


0.128659(1) 


2x16x16x256 


419200 


0.128658(1) 



Table 1: Critical hopping parameter k c for various lattice sizes. 

Lattice: 2x16x16x128 
/c =0.1286585 



5000 




Figure 2: Distribution of the order parameter s\ og at the transition point. 



in the case of larger volumes (L xy = 12, 16). For larger volumes, the clear two-peak structure 
with the equal-height criterion for the si og -distribution has been used to determine the transition 
point (cf. figure 0). In the case of smaller lattices, k c has been provided by the peak in the 
susceptibility. 

For the determination of the smallest mass in the Higgs channel the operators L^ XIJi and p x 
have been used at the transition point. The correlation functions have been determined in the 
^-direction with a relative error of O(10~ 2 ). In the case of smaller lattices {L xy = 4, 8), mass fits 
with one exponential can be used, and the determination of the splitting mass does not cause 
any problem. However, for L xy = 12, 16 the analysis is more complicated. By use of a one-mass 
fit we have determined the masses for different z-slice intervals with length of 10 lattice points. 
Let rrii denote the mass, fitted on the ^-slice interval: i < z < i + 10. Since the role of the larger 
masses becomes less and less important for larger distances one obtains the following condition: 

if i > j ==^ rrii < mj- (4) 

For the real data, decreases for a while; then, due to the noise, it starts to increase. We have 
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L X y 


4 8 12 16 


Eq from L v 
Eq from p 2 


0.147(3) 0.066(4) 0.034(1) 0.0040(6) 
0.149(4) 0.076(4) 0.033(6) 0.0069(1) 



Table 2: The measured energy splitting for different lattice extensions. 



chosen our fit interval in such a way, that monotony holds. Typical results are 5-25 and 5-55 
for L z = 128 and L z = 256, respectively. For these intervals we have performed two-mass fits. 
Since there are several local minima in the x 2- fit, we take those results, which appear for both 
z-extensions L z = 128 and L z = 256. The statistical errors are estimated by a jackknife analysis 
with 256 independent subsamples. The results are listed in table 2. The origin of the difference 
between the E values for L xy = 16 comes from the following fact. In the studied region the 
decay of the L v correlation function is dominated by the smallest mass, however, in case of p 2 
the decay is still dominated by the next excitation (E\ ~ 0.1). Therefore, it is better to consider 
Eq from p 2 as an upper bound. Figure |3] contains — \og(E )/L t as a function of I? . The slope 
of the fitted line gives the interface tension 

L = 0.053(5). (5) 

c 

The error contains the statistical error and an estimate of the systematics. The latter has been 
obtained from the difference between the results of the two different operators L v ^ and p\. The 
simulations have been performed on the Quadrics-APE machine of DESY, and as a consistency 
check, some of the calculations have been repeated with 64-bit floating point arithmetics. 
The prediction of the perturbative approach || is 

^ = 0.060(6). (6) 

c perturbative 

The error displayed here comes from the uncertainties in the lattice determination of the renor- 
malized masses and coupling ||. 



4 Two-coupling method and bootstrap analysis 

This method is described in ref. ||, therefore we just sketch the basic ideas. 

The lattice with z-extension L z much larger than the others, is divided into two halves [19 
In the part with lower ^-coordinates we choose the hopping parameter to be K\ < k c and in 
the other one to be K2 > K c . With a suitable initialization a mixed state with two interfaces is 
obtained. Let L^>(ki, k 2 ) denote the expectation value of L v in the part with hopping parameter 
«j. One can estimate the interface tension by 

(7~- lim lim ((K2-K1) I™ L z [.^(/Ci, «*) - ^L 2) («i, j ■ (7) 

Z K.2\K C Kl/Kc I. L z ^OO L r J J 



6 




Ki k>2 sweeps 




0.12845 0.12885 10000 
0.12849 0.12881 10000 
0.12853 0.12877 10000 
0.12855 0.12875 12000 
0.12857 0.12873 14000 


1.0607(30) 7.9042(91) 6.8435(83) 
1.0614(32) 7.007(10) 5.9461(92) 
1.0587(41) 6.058(12) 4.999(10) 
1.0684(38) 5.581(12) 4.513(11) 
1.0678(42) 5.035(12) 3.967(11) 



Table 3: Results for W, and AL^ on a 2 x 24 2 x 200 lattice. 



With the ansatz 

K 2 ) = — + bi + di ■ («j - K c ) + 0(Ki - K c f (8) 

a finite volume estimator of the interface tension is obtained from 

b = L z ■ (ci + c 2 ) . (9) 

As eq. (0) illustrates, a depends only on the difference AL V = [Ly — L^]. 

In table | the results for L$\ and AL^ on a 2 x 24 2 x 200 lattice are listed. We estimated 
the integrated autocorrelation time from the quotient of true and naive statistical errors. For 
and L)p the autocorrelation time is always found to be smaller than 50 sweeps. The errors, 
quoted in the table, are obtained by binning. 

The statistical errors of AL^ are always smaller than those of L@\ This is due to a kind of 
shift of the interfaces during the simulation; for instance, in the case of an interface moving from 
the region with k 2 into the one with Ki, both LM w ill increase, but the difference between them is 
not that much affected. Thus, there is a strong correlation between the two L^- values obtained 
on the same configuration. This has to be taken into account when estimating statistical errors 
for the interface tension, and is achieved by a bootstrap analysis. 
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lattice: 2x24x24x200, ,3 = 8.0, A = 0.0Q03 
_ 080 I 



1 .075 ; 
1 .070 ; 
1 .065 - 



1 .060 - 



1 .055 



1 .050 



0.12840 0.12845 0.12850 0.12855 0.12860 
k<k =0.1286565 



lattice: 2x24x24x200, (3 = 8.0, A = 0.0003 

9 . t V'd 



8.0 



7.0 



6.0 t 



5.0 



4.C 



0.12870 0.12875 0.12880 0.12885 0.12890 
k>k =0.1 286565 



Figure 4: Least x 2 ~fits for the determination of o from the two-coupling method. The x 2 ~ s are 
Xi = 1-88 and x\ — 1-50, respectively. 



The bootstrap procedure has been discussed in and used e.g. in We shortly discuss 
our analysis. Consider a set of N independent measurements for some primary quantities from 
a single Monte Carlo run. To form a bootstrap sample from the original N measurements we 
take N objects of them with repetition, thus some may be chosen a couple of times. The number 
of bootstrap samples is huge, N N , therefore one randomly generates representatives. Since 
correlations between the different primary quantities should be considered, these quantities have 
to emerge from the same Monte Carlo configurations for a given bootstrap sample. One computes 
in each sample a given secondary quantity a from the sample averages of the primary quantities 
(e.g. mass fits from correlation functions). An estimate of the statistical error is given by the 
dispersion of the a-values, thus by the two values left and right to the median containing 68.3% 
(34.15% in both directions) of the distribution. 

In order to eliminate the influence of autocorrelation we have combined 200 measurements to 
bins for each K-pair. (E.g. for the first row of table [3] the bin averages provide N = 50 independent 
values for the quantities and .) Our data are not obtained in a single Monte Carlo run 
but in 5 different runs for different K-pairs. For each /t-pair we generated a bootstrap sample 
and for the obtained samples we determined L)p and L$ as a function of k. Using the form 
of eq. (]|) we performed for these functions two x 2 -Ats. This procedure has been repeated 10000 
times. The statistical error has been determined from the distribution of the o. 

Figure |4] contains the combined data points with the corresponding fits. This yields 

-% = 0.065(9 + 1). (10) 

c 

The first number in the parenthesis is the statistical error, the second one stands for the uncer- 
tainty in the critical hopping parameter. For the individual fits we have c% = 0.00000277(290 + 5) 
and C2 = 0.0000378(61 + 7) with x\ — 1-88 an d x\ — 1-50, respectively. The error on (L 2 T~ 3 • c-i) 
is even larger than that of <r/T c 3 . This is an indication that the correlations have been taken into 
account. In order to check our error estimate based on the bootstrap procedure, we have also 
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performed a jackknife error analysis with 80 independent jackknife samples. The results of the 
two different methods are in complete agreement. 

In order to check whether the measured curvature of L^'(ki, k 2 ) as a function of is really 
described by eq. (H), we have also tried a 3-parameter quadratic ansatz. The results are x\ — 1-85 
and \\ = 4.15. The value of \\ shows that the data cannot be well described by a quadratic 
fit. Due to the large error bars in the symmetric phase (cf. figure |]) x\ is n °t very sensitive to 
the proposed fit-function. A 4-parameter fit with a quadratic term in addition to eq. (H) results 
in xi — 1-85 and x\ = 1-49. Since the x 2_s are nearly unchanged, the fourth parameter does 
not seem to represent a new degree of freedom and the value for <t/T c 3 = 0.063 does not deviate 
significantly from that of eq. (|lT)|) . 

It is also possible to determine a by using the AL^-values of table 3 together with the rough 
estimate k c = 0.12865 for the transition point, which has been obtained by the two-coupling 
method. One gets a/T 3 = 0.075(11) with x 2 = 0.77. In this case the error comes from 1000 
normally distributed random data and the uncertainty in k c has not been considered. With 
individual fits for both Z,W as previously described, the result is <r/T c 3 = 0.075(11) for k c = 
0.12865 with x\ = 1-89 and x\ = 1-54. From this agreement we conclude that the correlations 
between the L^-values are treated similarly in both methods. However, the difference between 
this result and that of eq. (|I0|) emphasizes the importance of a precise /^-determination. 



5 Conclusions 

We have determined the interface tension of the SU(2)-Higgs model with two different methods. 
Within the error bars the results of eqns. @ and fllPf ) agree with each other and with the 
prediction of the perturbation theory. When comparing the errors, it should be kept in mind 
that the product of the lattice volume and the number of sweeps (thus, CPU-time) is an order 
of magnitude smaller in the case of the two-coupling method. 



Together with the results of ref. ||, in which Binder's method |22] has been used with mul- 
ticanonical simulation, in total three different approaches for the measurement of the interface 
tension were studied. No significant difference between the results of these methods has been 
observed. The two-coupling method turned out to be the most economic way to estimate the 
interface tension. 
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